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Abstract. We investigate the statistics of the mean magnetisation, of its large deviations 
and persistent large deviations in simple coarsening systems. We consider more specif- 
ically the case of the diffusion equation, of the Ising chain at zero temperature and of 
the two dimensional voter model. For the diffusion equation, at large times, the mean 
magnetisation has a limit law, which is studied analytically using the independent interval 
approximation. The probability of persistent large deviations, defined as the probability 
that the mean magnetisation was, for all previous times, greater than some level x, decays 
O I algebraically at large times, with an exponent 9{x) continuously varying with x. When 

a: = 1, 6'(1) is the usual persistence exponent. Similar behaviour is found for the Glauber- 

CN ! Ising chain at zero temperature. For the two dimensional Voter model, large deviations of 

QQ . the mean magnetisation are algebraic, while persistent large deviations seem to behave as 

t^^ ■ the usual persistence probability. 
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1 Introduction 

Up to now most studies of persistence in simple nonequilibrium systems have focussed 
on the behaviour of the persistence probabihty at large times and on the computation of 
the related persistence exponent [1-28] . The aim of this paper is to broaden the scope of 
these former studies by investigating the statistics of more general persistent events. We 
will see that consideration of these events leads in particular to the introduction of new 
nontrivial exponents. 

A simple definition of persistence may be given as follows. Let a time-dependent random 
variable a{t) take only two values ±1, with some dynamical rule. Think for instance of 
a{t) as being the spin at a particular site in a dynamical Ising model. The persistence 
of this random variable up to time t corresponds to the most extreme situation where 
it never changed sign. In other terms the spin spent all its time in only one of the two 
possible phases. Note that, by its very definition, this event is non local in time. The 
probability of this event, or persistence probability, for most of the systems mentioned 
above, decreases algebraically in time, with nontrivial exponents. The surprise of finding 
new nontrivial exponents in the dynamics of nonequilibrium systems motivated to a large 
extent the interest for the subject. 

In this paper we investigate the statistics of large deviations and of persistent large 
deviations for simple coarsening systems. Both are natural generalisations of the concept of 
persistence. We apply this study to the case of the diffusion equation, the one dimensional 
Glauber-Ising chain at zero temperature, and the two dimensional voter model. 

Let us define the mean 'magnetisation' at time t of the random process cr(t), or 'spin' 

for short, as 

1 '■' 



M{t) = - / dua{u), (1.1) 

t Jo 

which is such that —1 < M{t) < 1. The quantities considered in this work are the following. 
- We first define the distribution of the mean magnetisation by 

P{t,x) = P{M{t)>x). (1.2) 

This quantity measures the chance for the mean magnetisation to deviate from its average.^ 
Taking a; > and t — > oo such that x is much larger than the width of the probability 
density function of M, defines the regime of large deviations. P(t, x) is then referred to as 
the probability of large deviations. 



^ Hereafter we will only consider cases with zero average magnetisation, i.e. such that the 
average {o'{t)) of the spin (and therefore of M{t)) over histories is zero. 



- We then define 

R{t, x) = P{M{u) >xyu<t). (1.3) 

This quantity will be hereafter referred to as the probability of persistent large deviations. 
Persistence, as defined above, corresponds to the largest deviation such that M{t) = 1 
(assuming for instance that a = 1 initially). The persistence probability therefore reads 

R{t) = P{a{u) = 1, Vw < t) = P{t, 1) = R{t, 1), (1.4) 

thus appearing as a limiting case of the two previous probabilities. 

We find, for the models considered in this work, the following results in the long time 
regime. 

(i) For the diffusion equation, M{t) has a limit law, i.e., P{t,x) converges to a limit 
distribution Poo{x) when t -^ oo. Using the independent interval approximation, the 
moments of this limit distribution are computed analytically. Their behaviour at high 
orders, or the singular behaviour of Poo{x) for x -^ 1, are related to the persistence 
exponent 9. 

The probability R{t^x) is found numerically to behave as t~^^^\ with an exponent 9{x) 
varying continuously from for x = —1, to 9, the usual persistence exponent, for x = 1. 
(Section 3.) 

(ii) For the ID Ising model at zero temperature, similar behaviour is found, namely 

P(t, x) -^ P^{x) and R{t, x) ~ t'^^^). (Section 4.) 

(iii) For the 2D voter model, numerical simulations suggest that, in the regime of large 
deviations, P(t, x) behaves as t~^^^\ with an exponent continuously varying with x, and di- 
verging when X ^ 1. They also seem to indicate that i?(t, x) behaves as exp[— J(x) (Int)^]. 
(Section 5.) 

We devote the next section to further considerations on large deviations and persistent 
large deviations. A general discussion and generalisations will be given in section 6. 

2 Large deviations and persistent large deviations 

Let us comment on the definitions (1.2, 1.3). 

First it is obvious that large deviations reflect a persistence property of the process. 
Think for instance of an event such that M{t) takes a value very close to 1, corresponding 
to a very large deviation. This is even more true of a persistent large deviation which is a 
more constrained event. 



Let us define the occupation time of tfie pfiases (+) or (— ), i.e. tfie time spent in tlie 
0" = ± pfiase, by 

In otfier words, tfie mean magnetisation M{t) of a generic spin gives a measure of tlie 
fraction of time tliat tfiis spin spent in one of the two phases. Therefore, while a large 
deviation only requires that a generic spin was most of the time in the same phase, a 
persistent large deviation constrains the spin to fulfil this condition at all previous times. 
Finally persistence corresponds for the spin to staying always in the same phase, hence 
R{t) = P{Tj^ = t). (Assuming that initially a = —1 would lead to the definition R{t) = 
PiT-=t).) 

Let us now illustrate the previous definitions on the very simple example of a symmetric 
random walk on a one-dimensional lattice. We denote by a{t) the step done by the walker 
at the discrete time t, where a = ±1 with probability 1/2. Starting from the origin at 
time 0, the position of the walker at time t is given by Xli^(^)- Its average position is 
equal to 0. The quantity M{t) introduced above represents the mean speed of the walker. 

The law of the position of the walker is well known. At large times it is normally 
distributed around its mean, with a variance proportional to t. As a consequence, the 
density of M is peaked around x = 0, with a variance decreasing as 1/t. The probability 
of a large deviation, giving a measure of the chance for the walker to reach a position far 
away from the origin, is exponentially small, and is given by (see Appendix A) 

P(t,x)~e-*^^") (x>0,t>l), (2.2) 

where I{x) = (1/2) [(1 + x) ln(l + x) + {1 — x) ln(l — x)] is an entropy function. In other 
words, the law of large numbers holds, the mean M{t) converging to its average (M) = 0, 
when t —^ oo. The limit law of M is a delta peak centered at x = 0, and P(t, x) -^ 
Poo{x) = H{—x), where H{x) is the Heaviside function. 

The persistent large deviation M{u) > Xj\/u < t corresponds to a situation where the 
walker always had a mean speed larger than x, i.e., stayed to the right of the position xt, 
between and t. If x > 0, R{t, x) behaves at large times in a similar fashion as in eq.(2.2). 
If X < 0, R{t, x) has a limit Roo{x) when t ^ oo which is a decreasing function of x, with 
a discontinuity at every rational value of x [29]. In the marginal case a; = 0, it is easy to 
show that i?(t, x) ~ 1/v^, for t large. 

Finally persistence corresponds for the walker to always stepping in the same direction. 
The persistence probability is 

i?(t) =e-*'"2_ ^2.3) 



Note that /(I) = In 2. 

By analogy with the case of the random walk, we set, for the models studied in the 
present work, 

P(t, x) ~ e-"(*)-^(^) (x>0,t>l), (2.4) 

which defines a function a{t) characteristic of the temporal behaviour of large deviations, 
and an entropy function I{x), keeping the same notation as above (see Appendix A). In a 
similar fashion, setting 

i?(t,x)~e-''(*)^^") (t>l), (2.5) 

defines a function b{t) characteristic of the temporal behaviour of the persistent large 
deviations. 

The time dependence of a{t) and b(t) for the models studied in this work is summarised 
in Table 1 and will be discussed in the next sections. 

Let us mention some mathematical references relevant for this work. Occupation times 
have been studied for Markov processes [30], and for several infinite particle systems [31, 
32]. Large deviations for occupation times were studied in [33, 34, 35, 36]. These references 
will be commented upon in the course of the paper. We are not aware of previous references 
on persistent large deviations. 

3 The diffusion equation 

3.1 The independent interval approximation and the persistence exponent 

We first introduce definitions, and remind results, which will be needed in the next 
section. Consider the equation 

at(^(x,t) = V2(/)(x,t), (3.1) 

where (/)(x, 0) is gaussian, with zero mean. Here x denotes a point in d-dimensional space. 
The changes of sign, or zero crossings, of the field (/> at a given space point, occur at 
times ti, ^2, . . . ,tn, starting from some time origin, or in the variable r = Int, at times 

Define, for a given space point x, the process $ = (^/a/((/)^). This process is gaussian 
and stationary in the time variable r, i.e. its two time correlation function ($(ri)$(r2)) = 
[cosh(r2 — ri)/2]~'^/^ only depends on the difference | r2 — ri | [11, 12]. As a consequence, 
the autocorrelation of the process a = sign($) reads [11, 12] 

A{r) = (a(O)a(r)) = - sin"! ^\,,a/2 - (3-2) 
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(3.4) 

Let us denote by In = Tn — Tn-i the intervals between zero crossings in the r variable. 
Considering the intervals as independent reduces the zero crossing process to a renewal 
process, entirely described, in the stationary regime, by /(/), the probability density func- 
tion of intervals. For such a process the probability Pni^) of having exactly n zero crossings 
up to time r reads, in Laplace space, 

Pn{s) = ^r-\s), (n>0) 

^^ ^ (3.5) 

11-/ 

*'»<'•' = " 7(0 ■ 

Noticing that A{t) = Xl^o(~)"P^(''^)' ^^ads to the following relation between / and A 
[11, 12] 

^ ^ l-{l)s{l-sA)/2 ^ l-s^d/2Ia{s) 
l + {l)s{l-sA)/2 l + sy/dj2Id{sy 



with (/) = TTy/sTrf. 

The persistence probability R{t) is the probability that the field (/> at a given space point 
did not change sign up to time t. Equivalently it is the probability that a^r) did not flip up 
to time r, i.e. the probability of no zero crossing Pq{t). At large times it behaves as t~^, 
or as e~^'^ . As a consequence, at large /, /(/) ~ e~^\ the persistence exponent appearing 
as the rightmost pole of /(s), s = —9 [11, 12]. 

For example, in one dimension: 

h{s) = V2J dr e— ^^^^ J = V2{P{s + 1/4) + P{s + 3/4)) 

oo . ^ ^ X (3.7) 

" ^^ U + 1/4+P ^ s + 3/4 + pJ' 

with /i(0) = 2n. The function /9(x) is related to ^(x), the logarithmic derivative of the 
gamma function, by 

m = l[^{^)+H^)]- (3.8) 
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In two dimensions: 



e"^" 



(3.9) 



,.0 «+l/2 + p' 



with 12(0) = TT. The largest zero of 1 + s^Jd/2Id{s) is found to be at s = —9, with 
e = .12032797884 . . ., for rf = 1 and = .186221071297 . . ., for rf = 2. 

3.2 Statistics of the mean magnetisation 

Our concern, in this section, is the determination of the distribution of the random 
variable M{t), at large times. We denote by t, or by r in the logarithmic scale, the 
observation time and by A the 'backward recurrence time', i.e. the length of time measured 
backwards from r to the last crossing event before r: X = t — Tn- The probability 
distribution of A in Laplace space reads, in the stationary regime,^ 

i(s) = mf). (3.10) 

We have, assuming that a{t) = 1, 

M(t) = ^ / dua{u) = Ut - tn - {tn - tn-l) + ■ ■ ■) =1-2^, (3.11) 

t Jo t 



where 



e = - - ^^^ + ■ ■ ■ = e-^(l - e-'" + e-^"-'"-i -...)= e'^X,. (3.12) 



Assuming that a{t) = —1 leads to M{t) = 2^ — 1. Note that $, = T^/t, according to the sign 
of (^{t), i.e. ^ is the fraction of time spent in the 'wrong' phase (cf eq. (2.1)). The random 
variable Xn = 1 — e~^" + e~^"~^"-i — ■ ■ ■ obeys the recursion relation Xn = 1 — e~^"X„_i. 
It is therefore recognised as a Kesten variable [37, 38, 39, 40]. 
To summarize, in the limit t ^ oo, the three equations 

M = ±(l-20, (3.13a) 

i^e-^X, (3.13b) 

X = l-e-^X, (3.13 c) 



^ Let us note that the age of the system considered in [21] is just equal to t — tn- It is 
related to the scaling variable tn/t = e~^, the distribution of which is known in the case 
considered here. 



contain the relevant information for the determination of fM{x) = — dPoo(x)/dx, the 
distribution of M. Equation (3.13 c) should be understood as an equality in distribution. 

The determination of the probability density of the Kesten variable X (and hence of 
fiviix)) for any given distribution of intervals /(/) is known in general as a hard problem, 
and does not seem feasible in the present case. However, from this set of equations we are 
able to extract the following information. 

(i) We perform a local analysis of fmix), for x ^ 1. 

(ii) We compute the moments {M^) of fuix). 

(iii) We solve the set of equations (3.13 a, b, c) for the case of an exponential distribution 
of intervals, /(/) = 6e~^\ proportional to the tail of the true distribution given, in Laplace 
space, by eq. (3.6). 

Let us analyse the local behaviour of fM{x) in the persistence region x ^ 1 (a similar 
analysis would hold in the limit x -^ — !)• Then ^ -^ 0, i.e, A -^ oo, with X finite. 
These conditions define the persistence region, where s ~ —9. Therefore /(/) ~ ae~ 
and q{X) ~ ae~^^ / {1)6, where a is the residue of /(s) for s = —9. Consider the Mellin 
transform of the law of ^, (C^)- From (3.13b) one gets 

{n=q{s){X^), (3.14) 

where q{s) is the Laplace transform (3.10). In this regime, one has {^^) ~ b/{s + 9), with 
b given by 

^=^(X-»>. (3.15) 

By inversion of the Mellin transform eq. (3.14) one obtains the behaviour of the distribution 
of ^, hence that of M in the persistence region x ^ 1. One finds 

fM{x)^2-^-^b{l-xf-^ (x^l). (3.16) 

As a consequence, for large k one has 

{M^) ^ 2-%r{9) k-'^ . (3.17) 

Note that the determination of b requires that of (X~^), which is unknown. A numerical 
estimate of the amplitude b can nevertheless be given, as follows. Using the method given 
in Appendix B, we computed the numerical values of the first 50 moments of M in one and 
two dimensions, in the independent interval approximation. By extrapolating these results 
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Figure 1: Plot of k^ {M^) (/c = 2, 4, . . . ,50) for the 2D diffusion equation in the inde- 
pendent interval approximation, versus /c~^, showing the approach to the limit amplitude 
.8424. Arrow: value of the limit amplitude .8476 for the beta law (3.18). 

we find 2-^6 r(^) ^ .870 for d = 1, and 2-^6 r(^) ^ .8424 for d = 2. In figure 1 a plot of 
k~^ {M'^) versus l//c, for (i = 2, is given, showing the approach to the limit amplitude. 

We are naturally led to compare the distribution fuix) to a beta law on (—1, 1), with 
the same singular behaviour in the region a; ~ 1, 
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The even moments of this law are given by 

_ i?((fc + l)/2,g) 
^^ 5(1/2,6') ■ 

The odd moments are zero, by construction. At large orders, 

,er(i/2 + ^) 



lik^2' 



k- 



(A;>1). 



r(l/2) 
Defining the ratio of local amplitudes A by /m ~ Af^^^^, for x ~ 1, yields 



A — lim 



k — >CXD ^/j 



hB{e,e + i), 



(3.18) 



(3.19) 



(3.20) 
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(3.22) 
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using (3.17, 3.21). Hence using the numerical estimates of b given above, we find A ^ .982 
for (i = 1, and A ^ .9938 for d = 2. Table 2 gives the values of the first moments {M^), 
compared to the moments of the beta law (3.18). 

The resemblance of fM{x) to the beta law (3.18) is enhanced by the fact that the solution 
of (3.13 a, b, c) for an exponential distribution of intervals /(/) = 9e~^\ proportional to 
the tail of the true distribution /(/) given, in Laplace space, by eq. (3.6), is precisely given 
by (3.18) (see Appendix B). This demonstrates the dominance of the tail of /(/) for the 
determination of fAiix). 

We also computed the probability distribution of the mean magnetisation obtained by 
numerical integration of eq. (3.1), for d = 1. This distribution is also found to be very 

close to /Beta_ 

In summary, the mean magnetisation M{t) has a limit distribution P^{x) when t -^ oo. 
In other words there is no law of large numbers for the random process o"(t), and absence 
of ergodicity. This distribution is found to be extremely close to the beta distribution 
(3.18). As long as 6* < 1, the density fM{x) diverges for x -^ ±1. Therefore the most 
probable values of M are near —1 and 1, while the average (M) = 0. The function I{x) 
diverges when x ^ 1. This signals the crossover to the persistence regime. The function 
a{t) introduced in (2.4) is formally of order unity. By extension, we will still speak of large 
deviations when t — ^ oo and a; ~ 1. 

3. 3 Persistent large deviations of the mean magnetisation 

We performed numerical simulations of the diffusion equation eq.(3.1) in ID, for a system 
size equal to 10^, starting from a random initial condition. At large times one observes an 
algebraic decay of the probability of persistent large deviations R{t, x) of the form 

i?(t, x) ~ t-^(^) (-l<x<l), (3.23) 

which corresponds to the behaviour h{t) ~ Int for the function defined in (2.5). The 
exponent 9{x) is to be identified to J{x) defined in eq. (2.5). 

Figure 2 gives a plot of the probability of persistent large deviations R{t, x) for a; = — .8. 
The usual persistence probability R{t) is also plotted, for comparison. The third curve 
corresponds to R{t, x, y) defined at the end of section 6 (see the comment there). 

Figure 3 shows a plot of the exponent 6{x) for — 1 < a; < 1. The exponent varies 
continuously from 0, for x = —1, to the value of the usual persistence exponent ~ .121, 
for X = 1. (We recall that the value of ^ ~ .121 obtained by numerical integration of 
(3.1) is slightly larger than the value of the exponent obtained by the independent interval 
approximation [11, 12].) We will further comment on these results in section 6. 
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Figure 2: Persistence probability -R(t), probability of persistent large deviations i?(t, x), 
and R{t,x,y) (see section 6), for the ID diffusion equation. (System size=10^.) From 
bottom to top: R{t) (slope= -.121), i?(t, -.8) (slope= -.096), i2(t, -.8, -.8) (slope= 
— .065). In dashed: regression lines. 
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Figure 3: Exponent 9{x) for the ID diffusion equation. Arrow: the usual persistence 
exponent 9 ~ .121. 

Let us finally mention that similar results as those presented in figures 2, 3 are found in 
two dimensions. 

4 The Glauber-Ising chain 
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We studied the Ising chain at zero temperature with the foUowing dynamics [41]. On 
each site of the ID lattice, values of the spin a = ± are initially distributed randomly. At 
each time step a site is picked at random. The spin on this site takes the value of one of 
its neighbours, chosen at random. 

We performed numerical simulations on a system of size L = 10^. As for the case of the 
diffusion equation, P(t, x) has a limit distribution when t —>■ oo. This distribution is very 
close to a beta law corresponding to the persistence exponent 9 = 3/8. The analytical 
study of P(i, x) will be given elsewhere. In particular it is easy to understand why this 
distribution converges to a limit when t — > cxo. For instance (M^) = A(l), where A{s) is 
the Laplace transform of the autocorrelation function A{t) = (cr(O)a(r)) with respect to 
the logarithmic time r = Int (see eq. (B.7)) [44]. 

The probability of persistent large deviations R{t, x) decays algebraically, with an ex- 
ponent continuously varying with x (figures 4 and 5). For x = 1 the usual persistence 
exponent 9 = 3/8 is recovered. 
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Figure 4: Probability of persistent large deviations R{t, x) for the ID Ising model. (Sys- 
tem size=10^.) From bottom to top: R{t), R{t, .5), P(t,0), i?(t, -.5), R{t, -.8). 



5 The 2D voter model 

The voter model is defined as follows [42]. On each site of a (i— dimensional lattice, 
opinions of a voter or values of a spin a = 1, 2, . . . , g are initially distributed randomly. At 
each time step a site is picked at random. The voter on this site takes one of the opinions 
of its 2D neighbours, with equal probabilities. Hence the rules of the dynamics of the voter 
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Figure 5: Exponent 9{x) for the ID Ising model. Arrow: the usual persistence exponent 

e = 3/8. 

model are a simple generalisation of those of the ID Ising model at zero temperature. In 
particular the ID voter model is identical to the Glauber-Ising chain. 

Earlier references to P{t, x) or to the occupation time T_|_ for the voter model may be 
found in [32, 34, 35, 36]. The function a{t) appearing in the large deviation expression eq. 
(2.4) is related to the variance of the occupation time T+ [32, 34, 35], hence to the two 
time correlation function of the process [44], by (see Appendix A) 



a{t) 



t' 



1 



VarT. 



+ 



VarM 



(5.1) 



In one dimension, VarM = const., as mentioned in section 4, hence a(t) = 0{1). The 
convergence in distribution of T^/t in one dimension was shown in [32]. In two dimensions, 
VarM ~ 1/lnt [32, 34, 35, 36, 44], hence a(t) ~ Int. Therefore the rate at which the 
distribution of M gets peaked is very slow. It was conjectured in [34, 35] that this estimate 
is exact, i.e. that a{t) = Int, hence that, in (i = 2, P{t,x) should have algebraic decay 
when t -^ oo. For d> 2, a{t) is respectively equal to -\/t, t/lnt, t for (i = 3, 4 and (i > 4 
[32, 34, 35, 36, 44]. 

We performed numerical simulations of the 2D voter model, for system sizes up to 
(4000)^, with q = 2. These simulations suggest that P{t^x) behaves as t~^^^\ with an 
exponent continuously varying with x (x > 0), and to be identified to I{x) defined in eq. 
(2.4). We will present the numerical analysis of the scaling in [44]. 

The numerical results also seem to indicate that R{t,x) behaves as exp[— J(a:) (Int)^], 
reminiscent of the behaviour of the usual persistence probability [16, 43]. Hence the func- 

13 



tion b{t) introduced in (2.5) is equal to (Int)^. It is nevertheless rather hard to conclude, 
on the basis of our numerical simulations. 

We conjecture that b{t) is proportional to N(t), the average number of particles in the 
dual particle system (diffusing coagulation, or equivalently reaction diffusion A+A^A, 
with a local source). N{t) is equal respectively to Int, (Int)^, v^, t/lnt, t for (i = 1, 2, 3, 4 
and d > 4 [32, 34, 35, 36, 43]. These results have a clear intuitive interpretation. When 
the dimension of space increases, particles interact less strongly since they have more and 
more space to explore before meeting. As a consequence, they are less correlated and 
their average number increases. In high enough dimension, the reaction between particles 
becomes irrelevant, hence N{t) becomes proportional to t, reflecting the total independence 
of the particles. Note that above two dimensions, a(t) is equal to N{t) (and therefore to 
b{t)). 

The voter model therefore interpolates between the case of the Glauber-Ising chain seen 
in previous section, if d = 1, and the case of the random walk of section 2, if d > 4. 

6 Discussion and conclusion 

The most striking conclusions that may be drawn from this study are, from our point 

of view, 

(i) the existence of a limit distribution for P(t, x), the probability distribution of the mean 
magnetisation, for the diffusion equation in all dimensions and the ID Ising chain, 

(ii) the appearance of families of exponents in the temporal decay of R{t,x), the prob- 
ability of persistent large deviations, for the diffusion equation and the ID Ising chain. 
Finally there are indications that the large deviations P{t, x) for the 2D voter model are 
algebraic. This extends the scope of former studies on the persistence exponents found in 
coarsening systems. 

A number of comments are in order. 

The time dependence of a{t) and 6(t), defined in (2.4, 2.5), is summarised in Table 1. One 
observes that for the random walk, a{t) and b{t) are proportional to t, while these functions 
are slowly increasing with time for the coarsening systems considered in this work, or even 
constant. This may be interpreted as follows. In the case of the random walk, M{t) is 
given by a sum of t independent random variables. In the coarsening systems studied here, 
values of the spin at a fixed position at different times are strongly dependent random 
variables. Thus a{t) and b{t) measure in some sense the effective number of independent 
variables in the system. This is clear for a(t), from its very definition given in Appendix A, 
and for 6(t), at least for the voter model, from the discussion given at the end of section 5. 
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Thus the mechanism by which R{t, x) -and as a consequence R{t), the usual persistence 
probabihty- have algebraic decay at large times, for the diffusion equation and the ID 
Ising model, can be traced back to the logarithmic behaviour of b{t) defined in eq. (2.5). 
The same comment holds for P(t, x) and a{t) for the 2D Voter model. 

Unfortunately the exact computation of the exponents seems a difficult task, since this 
amounts to computing the 'entropy' functions I{x) or J{x). Already computing the usual 
persistence exponent, corresponding in the present framework to taking the limit x -^ 1, is 
in general difficult. Moreover, even for the simple random walk, the probability R{t, x) is a 
nontrivial mathematical object. At least for the diffusion equation, for which it is possible 
to get analytic results for P{t, x) in the independent interval approximation (see section 
3), one could hope of computing R{t,x). Let us note that some aspects of this work, for 
instance the interpretation of the exponents as entropy functions, are reminiscent of the 
multifractal formalism. 

One also observes that for coarsening systems, when x ^ 1, I{x) diverges, while J{x) 
converges to a constant. The divergence of I{x) signals the crossover from large deviations 
{P{tjx)) to persistence {R{t)). The convergence of J{x) shows that R{t^x) is a natural 
generalisation of the persistence probability -R(t), h{t) encoding the type of decay of the 
persistence probability, being algebraic or not. 

We can enhance the difference between R{t, x) and P(t, x) as follows. First define the 
new random variable 

(j(t, x) = sign(M(t) -x), (5.1) 

which is an indicator of whether the mean magnetisation M{t) at time t is above or below 
the level x. One has 

P(t..) = (l±J^). (5.2) 

Then R{t,x) = P{M{u) > x,\/u < t) is just the persistence probability of this random 
variable. Therefore 

;j((,,).p(,(„,,).l,V„<t) = (i±i!l!^i±^...i±^), (5,3) 

(taking a discrete set of intermediate times, then letting n -^ oo) which shows that R{t, x) 
is a highly non local function of time. 

One may generalise the present approach by progressively 'thinning' large deviations, 
and tracking rarer and rarer events. Let us define 

1 <■' 



M{t,x) = - dua{u,x), (5.4) 

t Jo 
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(—1 < M(t,x) < 1) and the corresponding probabilities P{t^x^y) = P{M(t^x) > y) and 
R{t, X, y) = P{M{u, x) > 2/, V-u < t). 

First consider the random walk. Take x = for simplicity. Then M(t, 0) is simply 
related to the fraction of time the walker spends on the right side of the origin. The limit 
distribution of this quantity (when t -^ oo) is given by the arc sine law [45]. 

We computed R{t, x, y) on the diffusion equation, and on the Ising chain. For example, 
figure 2 shows R{t, x, y), with x = y = —.8, for the ID diffusion equation. Again algebraic 
decay is observed. 

Let us point out that this progressive thinning of large deviations implies probing the 
system by events which are always more non local in time. This questions the possibility 
of the existence of an infinite number of exponents in temporal quantities measured on 
strongly interacting systems. 
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Appendix A 

This appendix provides an explanation of the large deviation expressions (2.2) and (2.4). 

Independent random variables 

Consider Y = tM = X]i^(^)' where the a are independent identically distributed 
random variables. The generating function of moments of Y is 

(e^^) = /y(s)=e*^(^), (A.l) 

where /y(s) is the Laplace transform of the density of Y, and K{s) is the generating 
function of cumulants of the random variable a. Inverting the Laplace transform yields 

friy) = [ ;^e-^+*^(^) = / ^e-*[^--^(^)], (A2) 

J 2m J 2m 

where x — y/t. For t — i> cxo we use the saddle point method to evaluate the integral. At 
the saddle point, K'{sc) = x. Defining 

I{x)=s,x-K{s,), (A.3) 

yields fviv) ~ Q-^-^i^) ^ Finally 

P(t,x) ~e-*^(^\ (A.4) 

which is (2.2). Note that I{x) is the Legendre transform of K{s). 

Let us apply this general formalism to the case of the random walker (see section 2). 
Then (e'^'^) = coshs, and K{s) = In cosh s. At the saddle point x — K'{sc) — tanhsc, 
hence 

»e=2ln^-^. (A.6) 



Noting that cosh Sc = (l-x^)"^/^, leads to /(x) = (l/2)[(l+x) ln(l+x) + (l-x) ln(l- 

Correlated random variables 
Consider the generating function of the cumulants of M, denoted by Cm 



X] 



i^M(«)=ln(e^^) = V^.^ (A.6) 



Assume that, when t — * oo, the c^ scale as [34] 

b 



[a{t)Y 
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J {t^ oo), (A.7) 



where a{t) diverges with t, and the b^ are constants. Hence a{t) ~ 1/ VarM. 
We now consider Y = a{t)M. Then 



^Me^^) = -i^KMMt))^5^^.- = K(.) (t^oo). (A. 

(^^sY^^^a{t)K{s)^ (A.9) 



a{t) ^ ' a{t) ^ ^" ^^ n\ 
Hence at large times 



from which one gets 

J 27ri J 27X1 

where x = y/a{t). Continuing as above, we obtain, at large times 

P(t,x) ~e-"W^(^\ (A.ll) 

which is (2.4). Again, I{x) is given by (A. 3) and K'{sc) = x, 

Here, the role of a{t) parallels that played by t for the former case of independent 
variables. This function can be therefore interpreted as the effective number of independent 
variables, in the case where the spins a at different times are correlated. 

Note that in the cases of the diffusion equation or the ID Ising-Glauber chain, a{t) = 
(9(1), because all cumulants become constant when t — > oo. 

Appendix B 

In this Appendix we first show how to compute the moments of the distribution of the 
mean magnetisation for the diffusion equation. We then solve the set of equations (3.13 a, 
b, c) for the case of an exponential distribution of intervals, /(/) = 9e~^\ proportional to 
the tail of the true distribution given by eq. (3.6). 

Using eqs. (3.13 a-c), the moments {M^) can be computed recursively as follows. The 
computation is done in three steps. 

(i) From (3.13 c) one computes the moments of X from those of e~', i.e. as functions of 
the coefficients /jt, recursively, where fk denotes /(s) for integer values of the argument. 
For instance 



l + /i 

{X-) = 



2. 1 - /l 



(l + /l)(l-/2) (B.l) 



.3, 1-2/1 + 2/2-/1/2 



(l + /l)(l-/2)(l + /3) 
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(ii) From (3.13 b), and using eq. (3.10), one has 

{e) = {e-"'){X'')=^^{X''). (B.2) 

(iii) By symmetry, only even moments of M are non zero. From (3.13 a) they are related 
to those of ^ by a binomial expansion. 
Finally one gets 

2N 1 2 1 - A 



(M^) = 1 - 
(M^) = 1 



(01 + /V 

8 (l-2A + 2/2-/i/2)(l-/3) (B.3) 

3(0 (i + A)(i_A)(i + /3) 



Replacing f{s) by its expression in terms of A{s) (see eq. (3.6)), permits to recast (B.3) 
into 

,j^4. ^ 1 _ (l-3ii + 4i2)(l-3i3) (B-4) 

l-2i2 
The first line of eq. (B.4) may be understood as follows. In the long time regime, using 
the logarithmic time r = Int, one has 

{M^{r)) = 2e-2" / drae"^ / ' drie"iA(r2 - n), (B.5) 

Jo Jo 

where A{t) is the autocorrelation function (3.2). Laplace transforming both sides of (B.5) 
gives 

f,re-'^MHr))^'-§±^. (B,6) 

hence, when t ^ oo, 

(M2)=i(l), (B.7) 

since the rightmost pole of the right hand side of (B.6) is at s = 0. The result (B.7) is 
generic for coarsening systems, whenever the autocorrelation function is scaling in the two 
time variables. In particular it holds for the Ising chain studied in section 4. 

We now show that the solution of eqs. (3.13 a, b, c) for /(/) = 9e~^'' is the beta law 
(3.18). Setting Z = e~^ in (3.13 c), the integral equation for the invariant distribution fx 
reads 

fx{x) = [ dz [ dyfziz) fx{y) d{x-l + zy), (B.8) 

Jo Jo 
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where fz, the probabihty density function of the variable Z, is known from that of the 
interval length /, /(/). For /(/) = 6'e~^' one has fzi^) = z^~^, which cast into eq. (B.8) 
leads to the solution 

fx{x) = B-\e + 1, e) x' (1 - xY-\ (B.9) 

where B{6 ^1,6) is the beta function. Then computing 

(r) = (e-^^)(X^), (B.IO) 

with (e"^") = q{s) =e/{s + e), and {X') = B-\e + 1,9) B{e + s + 1,9) leads to (e) = 
B-^{9, 9 + l)B{9 + s,9 + 1) hence to the law 

f^{x) = B-\9, 9+1) x^-^ (1 - x)^, (B.ll) 

for the random variable ^. Finally, for this choice of /(/), the solution of (3.13 a, b, c) is 
a beta law on (—1, 1) 

/^^*^(x) = B-\l/2, 9) (1 - xy-\ (B.12) 

which is eq. (3.18). 
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a{t) 


m 


I{x ^ 1) 


J{x -^ 1) 


Random walk 


t 


t 


In 2 


In 2 


J Diffusion equation 
\ ID Ising 


0{1) 


Int 


oo 


e 


2D Voter 


Int 


(lnt)2 


oo 


const. 



Table 1: Summary of results for the functions a(t), 6(t), /(x), J{x) (see eqs. (2.4, 2.5)). 





(M2) 


^2 


(M4) 


/"4 


(M6) 


/"6 


ID 


.7996 


.8060 


.7383 


.7462 


.7035 


.7119 


2D 


.7268 


.7286 


.6459 


.6482 


.6008 


.6032 



Table 2: Values of the moments (M^) for the diffusion equation, computed in the inde- 
pendent interval approximation, compared to the moments of the beta law (3.18). 



23 



